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Abstract. 

Gravastar models have recently been proposed as an alternative to black holes, 
mainly to avoid the problematic issues associated with event horizons and singularities. 
In this work, a wide variety of gravastar models within the context of nonlinear 
electrodynamics are constructed. Using the F representation, specific forms of 
Lagrangians are considered describing magnetic gravastars, which may be interpreted 
as self-gravitating magnetic monopoles with charge g. Using the dual P formulation 
of nonlinear electrodynamics, electric gravastar models are constructed by considering 
specific structural functions, and the characteristics and physical properties of the 
solutions are further explored. These interior nonlinear electrodynamic geometries are 
matched to an exterior Schwarzschild spacetime at a junction interface. 
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1. Introduction 

The generally well-accepted notion of a black hole is deeply rooted in the general 
relativist community, although a considerable number of particle and condensed matter 
physicists view the respective notion of event horizons with some suspicion and alarm 
[H El E]. In addition to this, one still encounters the existence of misconceptions 
and a certain ambiguity inherent to the Schwarzschild solution in the literature (see 
Ref. [Ij for a review). Despite the fact that the evidence of the existence of black 
holes in the Universe is very convincing, a certain amount of scepticism regarding 
the observational data is still encountered [5]. This scepticism has inspired new and 
fascinating ideas. In particular, models replacing the interior Schwarzschild solution 
with compact objects, and thus, doing away with the problems of the singularity at 
the origin and the event horizon, have been proposed to some extent in the literature. 
However, we do emphasize that the interior structure of realistic black holes has not 
been satisfactorily determined, and is still open to considerable debate. In this line of 
thought, an interesting alternative to black holes has recently been proposed, namely, 
the "gravastar" (^fratitational vacuum star) picture developed by Mazur and Mottola [2]. 
In this model, and in the related picture developed by Laughlin et al the quantum 
vacuum undergoes a phase transition at or near the location where the event horizon is 
expected to form. The Mazur-Mottola model is constituted by an onion-like structure 
with five layers, including two thin-shells, with surface stresses a± and V±, where a 
is the surface energy density and V is the surface tangential pressure. The interior of 
the solution is replaced by a segment of de Sitter space, which is then matched to a 
finite thickness shell of stiff matter with the equation of state p = p. The latter is 
then matched to an external Schwarzschild vacuum with p = p = 0. Related models, 
analyzed in a different context have also been considered by Dymnikova [6]. 

Although Mazur and Mottola argued for the thermodynamic stability of their 
configuration, the full dynamic stability against spherically symmetric perturbations 
was carried out in Ref. [7j, through a simplified construction of the Mazur-Mottola 
model. The configuration was simplified to a three-layer solution, i.e., a de Sitter 
interior solution was matched to a Schwarzschild exterior solution at a junction 
surface, comprising of a thin shell with surface stresses a and V. It was found that 
many equations of state exist that imply the dynamical stability of the gravastar 
configurations. The latter dynamical stability was generalized to an anti-de Sitter or de 
Sitter interior and a Schwarzschild- (anti)-de Sitter or Reissner-Nordstrom exterior [8]. 
In Ref. [9], by disregarding the presence of thin shells, a generalized class of gravastar 
models was analyzed. It was found, using these models, that gravastars cannot be 
perfect fluid spheres, and that they necessarily exhibit anisotropic pressures. Specific 
stable solutions, with respect to axial perturbations, possessing continuous pressures 
were also further explored [10]. In this context, one may also replace the interior de 
Sitter spacetime with a solution governed by the dark energy equation of state, u = p/p, 
where u < —1/3 [11]. Note that the particular case of u = —1 reduces to the de Sitter 
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solution, and for < —1, one ends up with phantom energy, an exotic cosmic fluid which 
violates the null energy condition. Now, the notion of dark energy is that of a spatially 
homogeneous cosmic fluid, however, inhomogeneities may arise through gravitational 
instabilities. Thus, these inhomogeneous solutions, denoted as dark energy stars [3l [TT]. 
may possibly originate from density fluctuations in the cosmological background. Note 
that now the pressure in the equation of state p = cup is a. radial pressure, and the 
transverse pressure may be determined from the fleld equations, as emphasized in Refs. 
[T2] . in a rather different context. 

In Ref. [13], motivated by low energy string theory, an alternative model was 
constructed by replacing the de Sitter regime with an interior solution governed by a 
Chaplygin gas equation of state, interpreted as a Born-Infeld phantom gravastar. In this 
context, the extension of the latter Born-Infeld gravastar variation to rather more general 
nonlinear electrodynamic models shall prove extremely interesting, and this is, in fact 
the analysis of the present paper. It is interesting to note that nonlinear electrodynamics 
has recently been revived, mainly due to the fact that these models appear as 
effective theories at different levels of string/M-theory, in particular, in Dp— branes and 
supersymmetric extensions, and non-Abelian generalizations (we refer the reader to Ref. 
[13] for a review), and have enjoyed a wide variety of applicability, namely, in cosmology 
[T5] . Traversable wormholes have also been explored in nonlinear electrodynamics 
[T6l fT7] . It was found that regular pure magnetic evolving wormholes exist [IE] . However, 
static spherically symmetric and stationary axisymmetric traversable wormholes have 
been ruled out, mainly due to the presence of event horizons, the non- violation of the null 
energy condition at the throat, and due to the imposition of the principle of flniteness, 
which states that a satisfactory theory should avoid physical quantities becoming inflnite 
[T7] . In fact, the Born-Infeld model was inspired mainly on the principle of flniteness, in 
order to remedy the fact that the standard picture of a point particle possesses an inflnite 
self-energy [18|. Later, Plebahski presented other examples of nonlinear electrodynamic 
Lagrangians [19], and showed that the Born-Infeld theory satisfles physically acceptable 
requirements. 

It was in nonlinear electrodynamics that the flrst exact regular black hole solutions 
to the Einstein fleld equation were found [20l EH l22]- In this context, magnetic black 
holes and monopoles [23], regular electrically charged black holes of a hybrid type 
containing a magnetically charged core [23], and electrically charged structures with 
a regular de Sitter center [25], have been found. In this work, we are interested in the 
construction of a wide variety of gravastar models coupled to nonlinear electrodynamics. 
Using the F representation, we consider speciflc forms of Lagrangians describing 
magnetic gravastars, which may be interpreted as self-gravitating magnetic monopoles 
with a charge g. Using the P formulation of nonlinear electrodynamics, it is easier 
to flnd electric solutions, and considering speciflc structural functions we further 
explore the characteristics and physical properties of the solutions. We match these 
interior nonlinear electrodynamic geometries to an exterior Schwarzschild spacetime, 
thus avoiding the problematic issues related to the singularities and event horizons. 
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This paper is organized in the following manner: In section [21 we outline the field 
equations, and analyze some characteristics of the solutions. In section [3l we consider the 
F representation of nonlinear electrodynamics, and explore specific magnetic solutions 
in some detail. In section IH we outline the dual P formalism and extensively analyze 
electric gravastar solutions. In section O we conclude. 



2. Field equations 



We consider the action of (3 + 1)— dimensional general relativity coupled to nonlinear 
electrodynamics given in the following form 

^ = / v^fl^ + ^(^)l d'x, (1) 

J LlOTT J 

where R is the Ricci scalar, and L{F) is a gauge-invariant electromagnetic Lagrangian. 
The latter depends on a single invariant F given hy F = F'^'^F^iy/A [19], where the 
antisymmetric tensor F^,^ = A^^^ — A^ y is the electromagnetic field. In Einstein-Maxwell 
theory, the Lagrangian takes the form L{F) = —F/Att. However, we consider more 
general choices for the electromagnetic Lagrangians. Note that the Lagrangian may 
also be constructed using a second invariant G ~ F^y*F^^ , where the asterisk * denotes 
the Hodge dual with respect to Qi^^y. However, we only consider F, as this provides 
solutions that are interesting enough. 

Varying the action with respect to the gravitational field provides the Einstein field 
equations G^y = SnT^^y, where G^^y is the Einstein tensor, and the stress-energy tensor, 
Tf,y, is given by 

T,y = g,yL{F)-F,^Fy''LF. (2) 

with Lp = dL/dF. 

The electromagnetic field equations are the following 

{F^'^Lp).^^ = 0, {*Fn.,, = 0. (3) 

The first equation may be obtained by varying the action with respect to the 
electromagnetic potential A^. The second relationship, in turn, may be obtained from 
the Bianchi identities. 

Consider a static and spherically symmetric spacetime, in curvature coordinates, 
given by the following line element 

where a and jS are functions of r. 

Taking into account the symmetries of the geometry, the non-zero compatible terms 
for the electromagnetic tensor are 

F,y = 2E{x^) 5l 5:) + 2B{x^) ^If 5t^ , (5) 

so that the only non-zero terms are Ftr = E{x^) and Fg^ = B{x^). 
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The Einstein tensor components, in an orthonormal reference frame, (with c = G = 
1) are given by 

p-2/3 

Gu = -r- (2/5V + e^^ - l) , (6) 
G,, = ^ {2a'r - e'^ + l) , (7) 

G^^ = GgQ = —— -(3' + a + ra" + r{aY - ra (3' , (8) 

where the prime denotes a derivative with respect to the radial coordinate, r. Now, it 
is a simple matter, using equation ([2]) in an orthonormal reference frame, to verify that 
T^i = —Tff. Thus, using the Einstein field equation, and taking into account equations 
(E]) and ([7j), we have + = 0, which provides the solution a = —[3 + G. The 
constant G may be absorbed by defining a new time coordinate, so that without a loss 
of generality, we may consider C = 0. 

For notational and computational convenience, we consider the metric fields in the 
following form 

2m(r 



g2a(r) ^ g-2/3(r) 



1 - 



(9) 



r 

so that the Einstein tensor components reduce to 

Gti = -Gff = — , G^^ = Ggg = — — . (10) 

The function m(r) may be considered the quasi-local mass, and is denoted as the mass 
function. The stress-energy tensor, in the orthonormal reference frame, assumes the 
diagonal form T^,> = dia.g{p, Pr,Pt,Pt), in which p(r) is the energy density, Pr{r) is the 
radial pressure, and Pt{r) is the lateral pressure measured in the orthogonal direction to 
the radial direction. Thus, the Einstein field equations finally take the form 

p(r) = - Pr{r) = -L- E'Lp = - ^ , (11) 

pt{r)=L- , . B^Lf = -- . 12 

We will consider lE"^ Lp\ < oo and {B"^ Lp/i^r"^ sin^ 9) \ < oo, as r — >■ 0, to ensure the 
regularity of the stress-energy tensor components. 

One may now define the factor a'{r) as the gravity profile, as it is related to the radial 
component of proper acceleration that an observer must maintain in order to remain 
at rest at constant r, 6, 0. Note that the radial component of proper acceleration is 
given by = (1 — 2m{r)/r) a'{r) (see Ref. [26] for details). Thus, the convention used 
is that a'{r) is positive for an inwardly gravitational attraction, and negative for an 
outward gravitational repulsion [26]. For gravastars an essential condition is that they 
have a repulsive nature, i.e., a' < 0, which for the present case amounts to imposing the 
following condition 
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We emphasize that what is required from a spherically symmetric solution of the Einstein 
equations to be a gravastar model, is the presence of a repulsive nature of the interior 
solution, which is characterized by the notion of the "gravity profile". The interior 
solution is then matched to an exterior Schwarzschild solution which is outlined further 
ahead. 

We also explore the energy conditions, in particular, the weak energy condition 
(WEC), which is defined as Tj^ciU^U" > where is a timelike vector. The fact that 
the stress energy tensor is diagonal will be helpful, so that the following three conditions 
are imposed 

P>0, p + Pr>0, P + Pt>0. (14) 

From the first inequality, we verify that m' > is imposed. Note that p+Pr = 0, so that 
the second inequahty is readily satisfied. Relatively to the third inequality, consider the 
factor 

P + Pt = -(E' + ^^.BA L, = ^, (15) 
V sm 6 / Svrr"^ 

where for simplicity, N{r) is defined as 

N{r) = 2m' — rm" . (16) 
Thus, in summary, to ensure the WEC, we simply need to impose that 

m' > and N{r) > (17) 

throughout the spacetime. 

In the present case, the term p + pt is equivalent to the anisotropy factor which 
is defined a.s A = pt — Pr- The latter is a measure of the pressure anisotropy of the 
fluid comprising the gravastar. A = corresponds to the particular case of an isotropic 
pressure gravastar, which in our case reduces to the trivial case of E = B = 0, as may 
be verified from equation ( fT5|) . Note that A/r represents a force due to the anisotropic 
nature of the stellar model, which is repulsive, i.e., being outward directed if pt > Pr, 
and attractive if pt < Pr- 

We consider a cut-off of the stress-energy tensor at a junction radius a, much in 
the spirit of the original Mazur-Mottola gravastar model [2]. For instance, consider for 
simplicity that the exterior solution is the Schwarzschild spacetime, given by 

ds^= -i^l-'^^dt^+(l-^y\r^ + r\de^ + sm^ed(l)^). (18) 

M may be interpreted as the gravastar's total mass. In this case the spacetimes given 
by the metrics (jl]) and f|T8|) are matched at a, and one has a thin shell surrounding 
the gravastar. Note that to avoid the presence of an event horizon, we need to impose 
a > 2M. Using the Darmois-Israel formalism, the surface stresses are given by 
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where the overdot denotes a derivative with respect to the proper time, r, and the prime 
represents, for the present derivative with respect to a. a is the surface energy 

density and V the surface pressure (see Refs. [27] for details). The static case is given 
by taking into account d = a = 0. The total mass of the gravastar, for the static case, 
with the junction interface oq, is given by 



where rUg is the surface mass of the thin shell, and is defined as nig = Ana'^a. 

It is also important to mention that for a self-gravitating object (solution to 
Einstein's equations) to be considered as a possible alternative to a black hole, which 
was one of the original ideas behind the gravastar model, the surface redshift should be 
able to reach values that are higher than that of ordinary objects. The surface redshift 
is defined by Z = AA/Ag = Xo/Xe — 1, where AA is the fractional change between 
the observed wavelength, Aq, and the emitted wavelength, Ag. Thus, according to our 
notation, the surface value takes the form 



Now, for a static perfect fluid sphere the surface redshift is not larger than Z = 2 |28j . 
However, for anisotropic spheres this value may be larger [29]. For the present case of 
nonlinear electrodynamic gravastars, the bounds on the surface redshift, in principle, 
place restrictions on the characteristic parameters of the nonlinear electrodynamic 
theory, which can be analyzed on a case by case, but we will not pursue this here. 

3. F representation of nonlinear electrodynamics 

A number of nonexistence theorems prohibiting regular electrically charged nonlinear 
electrodynamic structures have been proposed. These basically state that any L 
approaching the Maxwell weak field limit, F <^ 1, such that L ~ —F, prohibit 
electrically charged static spherically symmetric geometries with a regular center |23j . 
However, it has been argued that due to the fact that the energy density attains a 
maximum at the center, the validity of the Maxwell field limit at the center cannot 
be expected [25]. Another important feature worth mentioning, is that imposing the 
Maxwell limit at the center and at infinity, inevitably leads to a branching of L as a 
function of F. However, this issue should not bother us, as in this paper we consider 
matchings of interior nonlinear electrodynamic solutions to an exterior Schwarzschild 
spacetime, and thus the range considered generally implies a monotonic behavior of L 
as a function of F. As the Maxwellian limit is not imposed at the center, but rather at 
infinity, the nonlinear electrodynamic constructions outlined in this work do not possess 
the weak field limit in the specific range of interest, i.e., < r < a, so that we are at a 
certain freedom to consider rather general Lagrangians, in particular, those that do not 
obey the Maxwellian limit for F <^ 1. 




(21) 



z 



g-a(ao) _ I 



(22) 
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Now, taking into account the electromagnetic field equations, equations ([3D, from 
{*F''^).^ = 0, we obtain B' = and E^g = 0, so that we have E = E{r) and B = 3(6). 
From (F^*^ Lp).^ = 0, we deduce 

ELf = ^, B = qmsm9. (23) 

The invariant F = F^^'^F^^/A, included for self-completeness, takes the following form 



Note that from the imposition of the WEC, equation (JTSll reduces to 



from which we verify Li? < 0. Now, from equations ( l23i) and ( l25ll . we determine the 
electric field, F, given by 

1 



F(r) = ^ -N{r)±\lN\r)-[^P^] . (26) 

One verifies that independently of N{r), the electric field diverges at the center. Thus, 
taking into account the principle of finiteness, which states that a satisfactory theory 
should avoid physical quantities becoming infinite, the message that one may extract 
is that for the present geometry, both electric and magnetic fields cannot coexist. One 
should either consider a pure electric field or a pure magnetic field. We emphasize that 
the principle of finiteness, which is a basic requisite of nonlinear electrodynamic theory, 
is related to physically measured entities or physical quantities, independently of the 
particular basis, as introduced by Born and Infeld and later retaken by Plabanski, such 
as the electric and magnetic fields. The F is a construct of the formalism, and therefore 
not the physically measured entity. 

3.1. Pure electric field 

Considering a pure electric field, B = 0, and the non-trivial case of the negative sign in 
equation (|26l) . we have 

N(r) o A^^ 



Sirqe ' 2(87rg, 

From ELp = qe/f"^, the following relationship is deduced 

r _ 8^9e 

and the Lagrangian is provided by equation ([T2|) . i.e., 



(28) 



L = -^^. (29) 
Sir r ^ ^ 

Relatively to the WEC, one may consider a brief analysis by taking into account the 
nonlinear electrodynamic quantities. From the condition p > 0, we deduce L < 2FLp, 
and as noted above, p + = is readily verified. From the condition p + > 0, we 
have FLp > 0, and as F < 0, then we verify that Lp < 0. 
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3.2. Pure magnetic field 

For a pure magnetic field, E = 0, we have the following relevant factors 

B = q^sme, P=^A^ (30) 



and 



(31) 



i. lib , . 

L = -7ZZJ- (32) 



Using equation (fTTI) . the Lagrangian is given by 

1 m' 
An 

From the above equations, namely, from F and L, one verifies that specifying a 
convenient nonlinear electrodynamic Lagrangian, equation (l32l) may be integrated to 
provide the mass function, and thus specifying the geometry of the solution. We explore 
several cases in some detail in the next section. 

Analogously to the previous case, one may analyze the WEC by taking into account 
the nonlinear electrodynamic quantities. From the inequality p > 0, we have L < 0. 
From the condition p + > 0, we deduce —FLp > 0, and as F > 0, then we verify 
that < 0. 

3.2.1. Magnetic Dymnikova gravastar. In this section, we analyze the Lagrangian 
analogue of the structural function, representing a regular electrically charged structure, 
proposed by Dymnikova [25]. Consider the Lagrangian and its derivative given by 

F , 1 



4n (l + b,/2^f ' 47r (l + b,/2qlF) 

where 6 is a characteristic parameter of the nonlinear electrodynamic theory. Note that 
the stability of this solution has also been analyzed in Ref. [30] . In the weak field 
limit, -F -C 1, the Lagrangian assumes the Einstein-Maxwell form, L ~ — F/(47r) and 
Lp ~ -l/(47r). 

Using F = q^/{2r^), the above relationships take the following form 

^ " ~87c{bg^ + r^y ' " ~47r(6^2 + ^2)3 ' (34) 

where we consider the definition g = \qm\- 

Now, using equation fl52]) . one may deduce the mass function, which is given by 

m(r) = -—tAt- — ^ + arctan ( ) . (35) 
A{bg^ + r^) AVb \Vbg J 

The stress-energy tensor components are given by 

9^ 5'^(r^ — bg"^) 

P = ~P' = 8n{bg^ + r^y ^* " 87r{bg^ + r^)^ ' ^ ' 



(33) 
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which are regular throughout the spacetime. The WEC is satisfied as can be verified 
from the following relationships 



m'(r) 



g2^2 



2^4 



which are manifestly positive. 

The metric fields are given by 



N{r) 



2g^r 



J2a 



-2/3 



+ 



(6c/2 + ^2)3 ' 



arctan 



(37) 



Defining y = r/{Vbg), the above equation provides 



arctan(?/) — 



y 



(38) 

(39) 
1.825, 



Using df/dy = 0, we verify that / possesses a single minimum, i/m 
independently of the value of b. Now, f{b,ym) = has a single positive root at 
b ^ 0.1775. For b > b^, we have f{b,y) > 0. If 6 < br, then f{b,y) possesses two 
roots, yi and y2, which refiects the existence of event horizons. Note that f{b, y) > for 
y < yi and y > y2, which are the cases we are interested in. This analysis is represented 
in figure [TJ 



0.5- 



f(b,y) 



-0.5- 




Figure 1. Roots of with bi = 0.2775, 63 = 0.1775 and 63 = 0.1075, 

respectively. Note the presence of two positive roots, for b < br — 0.1775, reflecting 
the existence of event horizons. For b > br, we have / > 0. 



Consider the gravity profile written in the following form 



G{b,y)=gb'/'a'{b,y) 



- arctan(y) 



4?/' ~ i [arctan(y) - ^] } 



(40) 



which is depicted in figured For b > b^, we verify that a' < 0, i.e., the gravity profile 
is negative, refiecting a repulsive character of the geometry, for all values of y. For the 
case of b < br and for / > 0, the gravity profile is only negative in the range y < yi and 
y > y2- In the right plot of figure Ej we have considered the specific case of 6 = 0.0975. 

Therefore, one may construct a nonlinear electrodynamic magnetic gravastar, by 
matching the solution outlined above to an exterior Schwarzschild solution. Note, 
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Figure 2. Plots of the gravity profile G{b, y) — g b^^'^a'{b, y). In the left plot, we have 
considered b > br, and have shown that a' < for all values of In the right plot, we 
have the case of 6 < 6^ and for / > the gravity profile is only negative in the range 
y < yi and y > t/2- We have considered the specific case of 6 = 0.0975. The dashed 
curves represent the gravity profile and the solid line depicts /(6, y). 



however, that the magnetic permeabihty, /i ~ ^/Lp, diverges at the center, r = 0, 
as Lp ^ 0, thus, exhibiting a magnetically superconductive behavior. For b > br, we 
verify the absence of event horizons and the gravity profile is negative for all values of 
y, so that the matching occurs simply at a > 2M. If 6 < br, then more care needs 
to be taken due to the presence of two horizons, at ri = Vbgyi and r2 = Vbgy2, 
respectively. The gravity profile is negative for r < ri, so that the matching should be 
at ri > ao > 2M, in order to avoid the presence of event horizons. The surface stresses 
of the thin shell are given by 



- . 1 




9' 



+ at. 



+ 



2Vb 



9 CLO 

arctan 



ao 



Vb9. 



4(692+^2) 



4v^ 



3— arctan (^) 



ao 



+ arctan (^) 



(41) 



(42) 



The total mass of the magnetic monopole gravastar is provided by the following 
relationship 



M 



9'^ao 



9 



^{bg^ + al) AVb 



arctan 



.Vb9. 



\ 



1 + 



2 (bg^ + al) 



2Vb 



arctan 



ao 



.Vbg^ 



2an 



(43) 



3.2.2. Magnetic Bardeen gravastar. Bardeen seems to have been the first author to 
surprisingly produce a regular black hole model [31]. It is interesting to note that 
the Bardeen model has been reinterpreted as a magnetic solution to the Einstein 
field equation coupled to nonlinear electrodynamics [32], i.e., it corresponds to a self- 
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gravitating magnetic monopole. (Gravitational magnetic monopole stellar solutions 
have also been explored in Majumdar-Papapetrou systems [33]). 

Consider the Lagrangian and its derivative, Lp, written in the following form 

5/2 



2^™^ ^ . 15 {2qjFy/^ ^^^^ 



Note, however, that this Lagrangian does not assume the Maxwell form in the weak 
field limit, i.e., L ~ —PF^^^ for F <til. It is perhaps important to emphasize here that 
to be a nonlinear electrodynamic model, the Maxwellian limit has to be recovered in 
the weak field limit. The nonexistence theorems, proposed in Ref. [23], of electrically 
charged regular structures imposed the weak field limit as r 0. However, it was 
argued that as the energy density attains a maximum at the center, the Maxwellian 
limit cannot be expected at r = Note that it is only at r — oo that the weak 
field limit is recovered F 0, for the spacetimes considered in this work. Therefore, 
in the gravastar models that we construct, it is not necessary to regain the weak field 
limit in the specified range < r < a, as the interior nonlinear electrodynamic solution 
is matched at a junction interface a, to an exterior vacuum geometry. 
Using F = g^/(2r^), then the above relationships take the form 

3^1 ^ 15 gsr'^ 



Airs (r2 + g%f^ ' Sirs (r^ + glf^^ 

where we consider the definition gs = \qm\- 

Using equation (152]) . one may deduce the following mass function, given by 

m(r) = — ^ . (46) 

s{r' + glf 

The stress-energy tensor components are given by 

3^1 _ 3^|(3r^-2 

Airs (r2 + glf ' * Sirs (r^ + gl 

The WEC is satisfied, as may be verified from the following relationships 

^'(0 = , iV(r) = . (48) 

s (r2 + gl) ' s (r2 + g^) 

The metric fields given by /(r) = 1 — 2m{r)/r, may be rewritten as 

using the definition y = r / gs- Taking into account df /dy = 0, we verify that / possesses 
a single minimum, ym = v^, independently of the value of s. Now, f{s,yrn) = has 
a single positive root at Sr ~ 0.7698. For s > Sr, we have f{s,y) > 0. If s < Sr, then 
f{s,y) possesses two positive roots, yi and y2, corresponding to two event horizons. We 
verify f{s, y) > 0, for y < yi and y > 2/2, which are the cases that we are interested in. 



(45) 



p=-pr= . ./%,5/2 > p^= r ), , (47) 
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The gravity profile given by 



may be rewritten as 

G = gsa' 



r2 + (?|)[2(7Br2-s(r2 + ^;|)3/2] 

y {y' - 2) 



2^3/2 



(50) 



(51) 



(1 + 1/2) [2y^-s{l + y 

using the definition y = t/qb- The case of s > s^, and for which a' < 0, is plotted 
in figure O Note that the gravity profile takes negative values, refiecting a repulsive 
character, in the range < y < ym for s > Sr- Thus, one may match this interior 
solution at = gs Vm > clq > 2M. For the case of s < and for / > 0, the gravity 
profile is only negative for y < yi, being positive for y > y2. Note that the gravity 
profile has asymptotes precisely at the roots of /. Thus, one may match this solution 
a.t ri = gsyi > ao > 2M, to an exterior Schwarzschild spacetime. 




0.5 



-0.5 













y2 










-- G(s,y) 








f(s,y) 



Figure 3. Plots of the gravity profile G{s, y) — gs a' , for the specific case of s > (left 
plot), and s < Sr (right plot), respectively. For the first case, s > Sr, the gravity profile 
takes negative values, reflecting a repulsive character, in the range < y < ym = V2- 
For the second case, s < we have used s = 0.5. The gravity profile is given by the 
dashed curves, and is negative in the range y < yi, where / > 0. The dashed curves 
represent the gravity profile, and the solid line depicts /(s, y). 



The surface stresses of the thin shell are provided by 
1 



a 



V 



47ran 



2Af 
ao 



1 - 



ao 



_ gB"oK+4g|) 



2M 
ao 



2g_Bao 



57rao I - ^ 
The total mass of the Bardeen gravastar is given by 

+ m^(ao) 



(52) 
(53) 



M 



gal 



s (ao^ + g^) 



2gal 



\ s (ao^ + g 



2an 



(54) 
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4. Dual P formalism 

Despite of the fact that nonlinear electrodynamics may be represented in terms of 
a nonlinear electrodynamic field, F^,^, and its invariants, one may introduce a dual 
representation in terms of an auxiliary field P^^,. The introduction of P^jy, proved to be 
extremely useful in the derivation of exact solutions in general relativity, especially in 
the electric regime. The dual representation of nonlinear electrodynamics is obtained 
by a Legendre transformation given by 

H = 2FLf - L . (55) 

The structural function if is a function of a factor P, defined by P = P^uP'^'^ /A. In this 
representation, the theory is reformulated in terms of the structural function H by the 
following relationships 

Pi^iu = LpF^y , F^y = HpP^y , (56) 
L = 2PHp - H , LpHp = 1 , (57) 
where Hp = dH/dP. The invariant P is given by 

The stress-energy tensor in the dual P formalism is written as 

T^, = g^, {2PHp -H)- P^^P,'' Hp , (59) 
and in the orthonormal reference frame, provides the following components 

= -T,, = H-^^P^^Hp, (60) 
sm d 

Tii = T^e = -H-PlHp. (61) 

The electromagnetic field equations in the P dual form are the following 

P'^^^ = 0, (ifp*P'^'^).^ = 0. (62) 

We emphasize that the tensor F^^ = Hp P^^ is the physically relevant quantity. 
However, one may obtain electric solutions easier by considering the P dual formalism. 

4-1. Electric field 

Equation (IH^ implies that for the pure electric field, P = 0, together with equations 
f l27|) - fl29|) provide the following extremely useful and simplified relationship 

He = ^%, (63) 
where Hp is the H structural function for this case and the correspondent P invariant 

P- = -il- (64) 
Thus, using the dual P formalism, it is easier to find nonlinear electrodynamic solutions 
than in the F formalism, for the specific case of pure electric fields. We consider several 
solutions in this section. 
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4- 1.1. Ay on- Beato— Garcia gravastar. In this section, we are interested in constructing 
a specific gravastar geometry from a regular black hole solution coupled to nonlinear 
electrodynamics, found by Ayon-Beato— Garcia [21]. Although this solution is indeed 
regular, i.e., the metric, curvature invariants and the electric field are regular everywhere, 
one still verifies the presence of event horizons, and consequently the associated 
difficulties related to these null hypersurfaces. Consider the specific Ayon-Beato— Garcia 
structural function [21] given in the following form 

-P 

H = . — , , (65) 

471 cosh^ [s^-2qlP ) 

where s is an adimensional constant. Note that H assumes the Einstein-Maxwell form 
in the weak field limit, i..e, H —P/At^ for P <^ 1. We also have 

1 



Using equation ([& 

,.2 

H 



8tt cosh^ (s^-2qlP 

, H and Hp may be recast as 
1 



s tanh^ [s{^-2qlP -2 



Svrr^cosh^ 
where the definition q 



— ] tanh^ 



' Svrcosh'f^^ 
is used. 

The mass function is determined by using equation ( l63ll . and is given by 



(66) 



(67) 



2s 



tanh 



(68) 



With this solution at hand, the energy density and the principal pressures take the 
following form 



P 



Pt 



Pr 



q^ 



iirr^ 



tanh^ ^ 



87rr^ 



tanh^ f — 



tanh ( — 



Analysing the geometry of the solution, consider the metric fields 



^2/3{r) 



1-^ 

s r 



tanh ( — 



Defining y 



rj (sq), the above equation takes the form 



s'^y 



(69) 
(70) 

(71) 
(72) 



Using df/dy = 0, which may be re-expressed as [1 — tanh(l/?/)] [y — 1 — tanh(l/?/)] = 0, 
we verify that / possesses a single minimum, y^ ~ 1.564, independently of the value of 
s. Now, /(s, ym) = has a single positive root at s ~ 0.5277. However, for this case, we 
have f{s,y) > if s > s^. If s < Sr, then f{s,y) possesses two roots, yi and y2, which 
refiects the existence of event horizons. 
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We also need to verify the WEC, given by m' > and > 0, which taking into 
account the definition y = r/{sq), assume the following form 



m{s,y) 
N{s,y) 



1 

g2y2 




(73) 
(74) 



which are represented in figure HI Note that at the center y = r = 0, the factor N 
becomes zero, showing that Pr = Pt, as was to be expected. We have only considered 
the range 0.3 < s < 1, for representational convenience. Note that m' > is verified for 
all values of s and y. However, N is negative in the range < y < yN — 0.4842, for all 
values of s. Thus the WEC is violated in the range < ?/ < ?/7v — 0.4842. 



m'(s,y) 




Figure 4. Plots of the WEC, m'(r) and N{r), respectively. We have only plotted the 
range 0.3 < s < 1, for representational convenience, to' > is verified for all values of 
s and y. However, N is negative in the range < y < — 0.4842, for all values of s. 



For the gravity profile we deduce the following relationship 



a 



|l — tanh 


f) 


-(?) 


1 - 


tanh^ (f )■ 


} 


2r 


sr 

. q 


— 1 + tanh 


(?)] 





(75) 



which using the definition y = r/{sq), takes the form 
G = {2sq)a' 



|l — tanh 


(i) 


-(i) 


1 - 


tanh^ (i)' 


} 


y 


s'^y 


— 1 + tanh 


m 





(76) 



The case of s > s,., and for which a' < 0, is plotted in figure [51 Note that the gravity 
profile takes negative values, refiecting a repulsive character, in the range < y < ym 
for s > Sr- Thus, one may match this interior solution at = sqym > cto > 2M. 
For the case of s < Sr and for / > 0, the gravity profile is only negative for y < yi, 
being positive for y > y2. Note that the gravity profile has asymptotes precisely at the 
roots of /. Thus, one may match this solution at ri = sqyi > qq > 2M, to an exterior 
Schwarzschild spacetime. 
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Figure 5. Plots of the gravity profile G = {2sq)a' , for the specific case of s > (left 
plot), and s < Sr (right plot), respectively. For the latter we have used s = 0.4277. 
The gravity profile is given by the dashed curves. For the first case, s > Sr, the 
gravity profile takes negative values, reflecting a repulsive character, in the range 
< y < Um — 1.564. For the case of s < s^. and for / > 0, the gravity profile is 
only negative for y < yi. 



Then, with our solution for the mass function, alongside H and Hp, we can cast 
the following relevant functions for the system 



E{r) 



Q 



l_tanh2 



2 I87rr2 
Lp = -Stt 

^2 



1 — tanh 
1 — tanh 



tanh 

2 fsQ 



r 

2 (SQ 



sq 



— tanh I — 



sq 



r 

2 (sq 



2 - 

sq 
r 



tanh 
tanh 



tanh 



sq 
r 



- 1 



2FLp = 



STrr^ 



l_tanh2^ 



sq 



— tanh ( — 



sq 



(77) 
(78) 
(79) 
(80) 
(81) 



Given this, the above interior Ayon-Beato— Garcia spacetime may be matched to 
an exterior Schwarzschild solution, at the junction interface, Oq, with the surface stresses 
given by 

1 / / 



a — 



47ran 



2M 
ao 



- Wl 



1 — tanh 



1 




- tanh f ^) + f ^) 


1 - 


-t-''=(S)l} 


87rao 


,/l-^ 

V »o 1 




1 — tanh ^ 







(82) 

(83) 



The total mass of the gravastar is given by 
1 - tanh ( — ) + ms(ao) S W 1 - 



M 



.2s 



1 

SttQ 



tanh 



2an 



(84) 
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where ms(ao) is the surface mass of the thin shell, given by ms(ao) = AiraQa. 

Summarizing the above results, we have the following: (i) For s > Sr, the factor 
/ is positive for all values of y, proving the absence of event horizons. However, the 
gravity profile, a', is only negative in the range < y < i/m, reflecting a repulsive 
character of the geometry, which is essential for gravastar solutions. Finally, the WEC 
is satisfied only in the range y > yisf, for all values of s. Thus, in conclusion, one may 
match this interior Ayon-Beato— Garcia geometry to an exterior Schwarzschild solution, 
at Tm = sqym > ao > 2M, whilst the WEC is violated in the range < y < yN < y-m- 
{ii) For s < Sr, the factor / has two roots, yi and ?/2- f{s,y) is positive for y < yi and 
y > y2, while the gravity profile is only negative in the range y < yi. As for the previous 
case, the WEC is violated in the range y < y^ < yi, N < 0. Therefore, to construct a 
gravastar geometry, one needs to match the Ayon-Beato— Garcia solution to an exterior 
Schwarzschild spacetime at ri = sqyi > oq > 2M, with the consequent violation of the 
WEC in the range < y < yjy < yi. Another important characteristic that both cases 
exhibit is an electrically superconductive behavior, at the center, as Li? ^ oo as r 0. 

One may also consider other regular black hole solutions obtained by Ayon- 
Beato— Garcia, in particular the structural functions given in Refs. [201122] . however, we 
shall not endeavor in this analysis. The message that one may extract, is that one may, 
in principle, construct a wide variety of gravastar models in the context of nonlinear 
electrodynamics, by using the regular black hole solutions found by Ayon-Beato— Garcia 
[201 [211 [22]. 



4.1.2. Nonlinear modified Tolman-Matese- Whitman mass function. Consider the 
following structural function and its derivative given by 

^ ^ ^52P(1 + W52P) ^ (l + Say/^) 

Note, however, that this structural function does not assume the Einstein-Maxwell form 
in the weak field limit, i..e, H ^ (— P)^/^ for P <^ 1. However, as mentioned above, we 
are not preoccupied with regaining the weak field limit in the specified range < r < a, 
as the interior nonlinear electrodynamic solution is matched at a junction interface a, 
to an exterior vacuum geometry. 

Taking into account P = — gg/(2r^), H and Hp take the form 

72 ' 



H= /"I., , Hp = -^^^^^L^—^ , (86) 



a^(l + ^) fl + ^ 

\ aqj \ aq 

where the definition q = |ge| is used. The mass function may be integrated to yield 

Arcr"^ 

m{r) = . (87) 

« (1 + -J 

Note that this mass function is similar in form to the Tolman-Matese- Whitman function 
analyzed in detail in Ref. [11], which represents a monotonic decreasing energy density 
in the gravastar interior, and was used previously in the analysis of isotropic fluid spheres 
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by Matese and Whitman [31]. Thus, we denote the mass function given by equation 
fl87|l the nonhnear modified Tolman-Matese- Whitman mass function. 
The stress-energy tensor components take the following form 

q{3aq + r^) {3aq - r^) 

P = -P- = —( ' = -( j_ 2^3 ■ (^S) 

a [aq + r^j [aq + r^j 

The WEC is also satisfied as the following factors 

// X 47rr2g (3ag + r^) Svrr^g (Sag + r^) 

m (r) = 2 — 5 ^v) = 3 — ) (89) 

a {aq + r^) a {aq + r^) 

are manifestly positive. Note that at the center, r = 0, we have = 0, so that p^ = pt, 
as was to be expected. 

Analysing the geometry of the solution, consider 

^2.(0 ^ ^-2m ^ 1 Svrr^ 

Defining = Snq/a and = r'^/{aq), so that the above equation takes the form 

f{s,y) = l-^^. (91) 

1 + r 

Note the absence of event horizons for s < 1. If s > 1, then / possesses a positive root 
situated at i/r = 1/ v? — 1 . 

The gravity profile is given by 

a' = ^^"^'^ . (92) 

{a'^q + ar"^ — Srcr'^q) {aq + r^) 

which using the definitions = Sirq/a and y"^ = r'^/{aq), takes the following form 

s^y 

G = y/aqa=-- „, „, ,^ ^. (93) 

^ [l + (l-s2)y2](i + ^2) 

One readily verifies that a' < for s < 1. For s > 1, we also verify that a' < for y < yr- 
See figure [6] for the latter case, with the specific choice of s = 1.1. An asymptote of a' 
exists precisely at the root. 

The electric field and the F invariant are given by 

_ {5aq + r^) p _ { r^baq + r^) ] ' 

~ ^ I 2^3 ' 9 ^ I 2^3 ' ^ 

respectively. We also have the following relationships 

^ _ g^ (Sag - r^) ^ ga (ag + r^)^ 

(ag + r2)=^ ' ^ (Sag + r^) ' ^ ' 

In summary, one may match the interior solution to an exterior Schwarzschild 
spacetime, at a junction interface, oq > 2M. If s > 1, then the matching needs to be 
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Figure 6. Plot of the gravity profile G = ^Jaq a', for the specific case of s > s^. 
We have used s — 1.1. The gravity profile takes negative values, reflecting a repulsive 
character, in the range < y < yr- 



done at = y/ctqyr > > 2M. The surface stresses are given by 
1 



Ana, 



1 







2M 
ao 



Svran 
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ina, 



(i + S 

The total mass of the gravastar is given by 



1 _ M 

2M 



Svrao 



(96) 
(97) 



M 



Anal 



(1 + 5) 



+ ms{ao) 



a' 



iTcan 



a2 1 + ^ 

' aq 



2ao 



(9^ 



where ms(ao) is the surface mass of the thin shell, given by ms(ao) = Ana^a. 
4-2. Magnetic field 

Now for the E = case, equations fl30|) - fl32|) together with fl55|) provide the following 
relationships 

1 m" 



Hb = -. 



57r r 



(99) 



where Hb is the H structural function for this case and the correspondent P invariant 
is given by 

Ar2 



We can, in principle, find new solutions for the mass function m(r) by choosing a 
suitable form of the H = Hb structural function. Although this treatment is far from 
being trivial, we have considered the F formalism to find purely magnetic solutions. 
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5. Summary and duscussion 

Gravastar models have recently been proposed as an alternative to black holes, mainly 
to avoid the associated difficulties with event horizons and singularities. In this 
work, we have been interested in the construction of gravastar models coupled to 
nonlinear electrodynamics. Using the F representation, we have considered specific 
forms of Lagrangians describing magnetic gravastars, which may be interpreted as 
self-gravitating magnetic monopoles with a charge g. Using the P formulation of 
nonlinear electrodynamics, it is easier to find electric solutions, and considering specific 
structural functions we further explored the characteristics and physical properties of 
these solutions. We have matched these interior nonlinear electrodynamic geometries 
to an exterior Schwarzschild spacetime at a junction interface a, thus avoiding the 
problematic issues related to the singularities and event horizons. We emphasize that 
what is required from a spherically symmetric solution of the Einstein equations to be 
a gravastar model, is the presence of a repulsive nature of the interior solution, which 
is characterized by the notion of the "gravity profile" . It is important to point out that 
to be a nonlinear electrodynamic model, the Maxwellian limit, L ~ —F and Lp ~ — 1, 
has to be recovered in the weak field limit, F <^ 1. For the spacetimes considered 
in this work, it is only at r cxd that the weak field limit, F — > 0, is recovered. 
Therefore, in the gravastar models that we construct, it is not necessary to regain the 
weak field limit in the specified range < r < a of the interior nonlinear electrodynamic 
solution. Thus, we have used general Lagrangians and structural functions that are 
strongly non-Maxwellian in the weak field limit, namely the Bardeen and the Tolman- 
Matese- Whitman solutions. In conclusion, a rather wide variety of gravastar solutions 
may be constructed within the context of nonlinear electrodynamics. 
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